Saturation of magnetic films with spin-polarized current in presence of magnetic field 
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Influence of perpendicular magnetic field on the process of transversal saturation of ferromagnetic 
films with spin-polarized current is studied theoretically. It is shown that the saturation current 
Js is decreased (increased) in case of codirected (oppositely directed) magnetic field and current. 
There exists a critical current Jc > Js which provides "rigid" saturation - the saturated state is 
stable with respect to the transverse magnetic field of any amplitude and direction. Influence of the 
magnetic field on the vortex-antivortex crystals, which appear in pre-saturated regime, is studied 
numerically. All analytical results are verified using micromagnetic simulations. 

PACS numbers: 75.10.Hk, 75.40.Mg, 05.45.-a, 72.25.Ba, 85.75.-d 



I. INTRODUCTION 

The influence of spin-polarized current on planar mag- 
netic systems is of high applied and academic interest 
now. It is so mainly due to the possibility to handle the 
magnetization states of magnetic nanoparticles (nano- 
magnets) without using the external magnetic fields of 
complex space-time configurations. That provides new 
opportunities in construction of purely current controlled 
deviceJ^, e.g. magnetic di sk d rivers or Magnetic Random 
Access Memory (MRAM^. 

A convenient way to provide the infiuence of spin- 
polarized current on the magnetic film is to use a pillar 
structure which was firstly proposed in Ref.l^. The sim- 
plest pilar structure consists of two ferromagnetic layers 
(Polarizer and Sample) and nonmagnetic Spacer between 
them, see Fig. [T] When the electrical current passes 
through the Polarizer the conduction electrons become 
partially spin-polarized in direction which is determined 
by the Polarizer magnetization. Polarizer is usually made 
of a hard ferromagnetic material whose magnetization is 
kept fixed. Spacer, being very thin (few nanometers), 
does not change spin polarization of the current elec- 
trons but it prevents the exchange interaction between 
Polarizer and Sample. Thus the spin-polarized electrons 
transfer the spin-torque from Polarizer to the Sample 
what can result in dynamics of the Sample magnetiza- 
tion. The spin-torque influence can be described phe- 
nomcnologically by adding the Slonczewski-Berger term 
into Landau-Lifshitz equatiorPl^. 

Recently we studied influence of strong spin-current on 
the magnetic films'^. It was shown that the strong spin- 
polarized current can saturate magnetic film and value 
of the saturation current density J, increases with the 
film thickness increasing. We also demonstrated that in 
the pre-saturated regime a stable vortex-antiv ortex lat- 
tices (VAL) appear. As it was recently showri^SIIll the 
external magnetic field can drastically modify the mag- 
netic system dynamics induced by the spin-torque. The 
aim of this paper is to study the influence of perpendic- 
ular magnetic field on the process of the film saturation 




FIG. 1: The simplest pillar heterostructure consists of 
two ferromagnetic layers: Polarizer and Sample, and 
one nonmagnetic layer between them (Spacer). Black 
(larger) arrow shows the current direction, which flow 
through the heterostructure, and white (smaller) arrow 
shows the magnetization direction of the Polarizer. 



with spin-current. For this purpose, we modify devel- 
oped in Ref.^ linear theory of instability of the saturated 
state for the case of presence of magnetic field and uni- 
axial anisotropy. It enable us to obtain the dependence 
of saturation current Jg on the field amplitude. We also 
demonstrate that in linear approximation the actions of 
the perpendicular magnetic field and uniaxial anisotropy 
on the stability of saturated state are equivalent. Using 
micromagnetic simulations we study how the properties 
of the VAL, which appear it the pre-saturated regime, 
depend on the value of the applied field. 



II. THE THEORY OF SATURATED STATE 
STABILITY 

Let us consider a ferromagnetic film with thickness h 
and lateral size L ^ h. We use here a discrete model of 
the magnetic media considering a three-dimensional cu- 
bic lattice of magnetic moments Mi, with lattice spacing 
a h, where u = a{vx,Vy,Vz) is a three-dimensional 
index with v^iVy^v^ € Z (here and below all Greek in- 
dexes are three-dimensional and Latin indexes are two- 
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dimensional). We assume also that h is small enough 
to ensure the magnetization uniformity along the thick- 
ness. In this case one can use the two-dimensional dis- 
crete Landau-Lifshitz-Slonczewski equatiorP^^ 



dE rrin x [m^ x z] 

— >c 

drrin 1 + (m„ • z) 



(1) 



to describe the magnetization dynamics under the influ- 
ence of a spin-polarized current which flows perpendicu- 
larly to the magnetic plane, along the i-axis, see Fig. [l] 
It is also assumed that the current flow and its spin- 
polarization are of the same direction in Eq. ([T]). The 
two-dimensional index n = a{nx,ny) with n,j;,ny G Z 
numerates the normalized magnetic moments m„ — 
Mn/\Mn\ within the film plane. The overdot indicates 
derivative with respect to the rescaled time in units of 
{A'K^Ms)~^ : 7 is gyromagnetic ratio, Mg is the satura- 
tion magnetization, and £ — E / {AnM'^a^J^z) is dimen- 
sionless magnetic energy, where Ji^ h/a is the number 
of magnetic moments along the thickness. The normal- 
ized current density is presented by parameter x = r\J j Jq 
, where 77 is the degree of spin polarization, J is the cur- 
rent density, and Jq = ATTM'^\e\h/h with e being electron 
charge and h being Planck constant. 

The Eq. ([T]) is written for the case when the conduc- 
tance of the Sample is much lower than the conductance 
of the Spacer, what corresponds to high level of spin accu- 
mulation at the nonmagnetferromagnet interfaces. The 
mismatch between spacer and ferrom agnet resistances is 
traditionally described by A-parametep^ES. But as it was 
shown in the Ref.l^ parameter A is not included in the lin- 
earized problem and therefore it has no influence on the 
saturation process, that is why we do not include A into 
our model assuming A ^ 1. We also omitted damping 
in the equation of motion ([T]), because, as it was shown 
earliei'^, the spin-current provides an effective damping 
which is much larger than natural damping. 

We consider here a magnetic system, the total energy 
E = E'ox + Ed + Ez + i?an of which consists of four parts: 
exchange, dipole-dipole, Zeeman and magnetocrystalline 
anisotropy contributions. Exchange energy up to a con- 
stant has the form 



E,^ = -S'^l^Jo 



rrir 



(2) 



where h numerates the nearest neighbors within the film 
plane of n-th atom, § is value of spin of a ferromagnetic 
atom, and 3o > denotes the exchange integral between 
nearest atoms. 

The energy of dipole-dipole interaction is 



E,, 



E 



(3) 

where rxv = \ — v with A and v being the three- 
dimensional indexes. 



The Zeeman energy describes the interaction of mag- 
netic film with external perpendicular magnetic field 
B = Bz and it reads 



Ez = -BM.a^yi^^; 



(4) 



And finally we introduce the energy of uniaxial 
anisotropy, which axis is oriented perpendicularly to the 
film plane: 



K 



(5) 



where K is the anisotropy coefficient which can be posi- 
tive (easy-axis) as well as negative (easy-plane) value. 
By introducing the complex variable^ 



mz, + im. 



y 



one can write the Eq. ([T]) in form 



di^n dtp* 



where function 



(6) 



(7) 



(8) 



represents an action of the spin-polarized current. 

For the future analysis it is convenient to proceed to 
the wave- vector representation using the two-dimensional 
discrete Fourier transform 




ik-n 



with the orthogonality condition 

J2 e*(fe-fc')-» = yi^yA{k - k') 



(9a) 
(9b) 

(10) 



where J^xy — E'^ ja^ is the total number of atoms 
within the film plane, fe = {kx,ky) = ^{lx,ly) is two- 
dimensional discrete wave vector, lx,ly G Z, and A(A;) is 
the Kronecker delta. 

Applying ([9| to the equation ^ one obtains equation 
of motion in reciprocal space: 



d£. 

Wk 



(11) 



Since we are studying the stability of the saturated 



state we can linearize the equation of motion (11) in 
vicinity of the solution = 1 what is equivalent to 
\ipri\ = and jT/ifcl = 0. To obtain the energy functional 
£ in -representation we expand components of the 
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magnetization vector into series in the way similar to the 
representation in terms of the Bose operators^"^: 



V2 



m 



(12) 



Substituting ( [T2| ) into energy components p| , (|3| , Q , 
([5]) and applying the Fourier transform (|9|, ( |10[ ) one can 
write the dimensionless energy functional in form 



Here the exchange contribution read^ 



(13a) 



where I — ^JW^^iocjJ^nM^a) is the exchange length with 
c = 4 being the number of nearest neighbors within the 
film plane. The energy of dipole-dipole interaction has 
the fornP 



£d = E 



+ 



g{kK) 



g(kh) 



(/fc^ - iWf 



fc2 



^h^-k + c.c. 



(13c) 



O(l^fcl^), 



where 17(0;) = (a; + e 
the form 



^ ~ l)/a;. The Zeeman energy takes 
£z = /3 5]|4|', (13d) 



where /3 = B/(4:TrMs) is the dimensionless magnetic field 
in units of the saturation field. And anisotropy energy 
can be written as 



(13e) 



where a — K/(AttM'^) is the dimensionless anisotropy 
coefficient. 

Details of deriving the contributions |13b| and |13c| in 
the wave-vector space can be found in Appendix A of 
the Ref.l2l. 

The current- action function 'J in the wave space has 
the simple form current: 



(14) 



Substituting (13) and (14 1 into (11) one obtains the 



system of linear equation for the complex amplitudes ipk 




FIG. 2: (Color online) The infiuence of parameter h on 
the instability domains (filled regions) of the 
transversally saturated Permalloy = 5.3 nm) film 
with thickness h = 20 nm. The instability domains are 



(13b) determined by condition (19). Points at the maximums 



determine the saturation current Kg for the given b. 



and V'lfc: 
where b = a - 



'kH^-l + i^+b + ^^ 
g{hk) (fc^ +ifcS')2 



/fc2 



k^r - 1 + ' + 
g(hk) (k'^ ^iky f - 

Vk- 



(15) 



Looking for solutions of Eq. ( 15 ) in the form 



(16) 



where are time independent amplitudes, one obtains 
the following condition for the rate constant z{k) 



z{k) 



± k(k). 



(17) 



Here the function x{k) is given by 



kik) = v/(l - fc2£2 _ 5) (p^2 _^ ^(/^fc) -l + b). (18a) 



Accordingly to the Eq. ( |17[ ) one can conclude that if 
value of the function >c{k) is complex then the saturated 
state of the film is stable. If value of the function >c{k) is 
real then wc have two different cases: for strong currents 
when >c > 2k we have Rez(fc) < and therefore the 
stationary state of the system is the saturated state with 
rriz — 1. However, for smaller currents >f < 2k the insta- 
bility of the saturated state develops. Functions 2k{k) 
for different values of parameter b is shown in the Fig. [2] 
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FIG. 3: (Color online) Dependence of the saturation 
current on the magnetic field. Solid line corresponds 
to the analytical solution obtained from ([Is]) and results 
of micromagnetic simulations (see text) are shown by 
markers. Dashed lines correspond to asymptotic (20b I 
and dot-dashed lines show the asymptotic (20a). 



Dependencies 1 and 2 corresponds to different 
thicknesses: h — 20 nm and h ~ 10 nm respectively. 



One can see that k{k) is a nonmonotonic function, which 
reaches its maximum value >Cc at k — K: 

=0, maxk(k) = k(Ko). (18b) 

dfc k 

Thus the value k^ determines the lowest current at which 
the saturated state remains stable - the saturation cur- 
rent: 



(18c) 



The instability domains, which are determined by condi- 
tion 



and X < Xs 



(19) 



are shown in the Fig. [2] as filled regions. Maximum val- 
ues of the shown dependencies determine the saturation 
current for the given value of the 6-parameter. De- 
pendence Xg (6) is shown in the Fig. [s] One can see that 
for 6 > 1 the magnetic film is perpendicularly saturated 
without current. For example, this case can be realized 
for magnetically soft film [a = 0) when the external field 
exceed the saturation value /3 > 1. Analysis of (18 1 en- 



able one to obtain the following asymptotical behaviors 



h fl-b 



3/4 



h^\h\ 



6<1, 



-oo. 



(20a) 



(20b) 



Accordingly to ( 20b I the saturation current is a bounded 
above quantity: Xs < 1 for any values of parameters. 
It means that for currents x ^ 1 the perpendicularly 
saturated magnetic film remains stable for any values of 
magnetic field and uniaxial anisotropy constant. In other 
words, if the current >f ^ 1 is applied then magnetization 
reversal is not possible with perpendicular magnetic field 
of any (even infinitely large) amplitude. We call this 
phenomenon "rigid saturation" . 

The critical current Xc = 1 has the following dimen- 
sional form 



47rM2|e|/i 
hrj 



(21) 



Thus the current is determined only by material pa- 
rameters (saturation magnetization) and thickness of the 
film. For the case of permalloy film with thickness h — 20 
nm and rate of spin polarization rj = OA the expression 
^ results Jc = 70.6 x lO^^ A/m'^. 



Thus we determine the physical meaning of the value 
Jo ~ the minimal current density which provides the rigid 
saturation (for the case of full spin-polarization rj = 1). 

To verify our analytical results we used full-scale 
OOMMjEl micromagnetic simulations. This modelling 
were performed with material parameters of permalloy as 
follows: saturation magnetization Ms = 8.6 x 10^ A/m, 
exchange constant A = 1.3 x 10~^^ J/m. These values 
of parameters correspond to the exchange length £ = 5.3 
nm and saturation field of the infinite film Bg = 1.081 T. 
Since the external field and the anisotropy are included 
i nto p roblem in equivalent ways, see Eqs. ( 13d[ ), ([s]), ( [T5| , 
(18a I, in the simulations we restrict ourselves only with 
case of magnetic field, neglecting the anisotropy (a = 0). 
Not being able to simulate he film of infinite size we chose 
two nanodisks with diameter D = 350 nm and thick- 
nesses ft, = 20 nm and h = 10 nm and mesh size was 
3 X 3 X ft nm. In the absence of the magnetic field and cur- 
rent the ground magnetic state of nanodisks of the men- 
tioned sizes is vortex distribution of magnetization, see 
Fig. [T] To these nanodisks we simultaneously apply the 



external magnetic field of form B{t) = Bq{1- 



-t/At 



) and 



spin polarized current J{t) — tAJ/At, where At = 1 ns, 
A J = 10" A/m^ with rate of spin polarization rj = 0.4. 
Gradual increase of the field and current allow us to avoid 
an intense magnon dynamics. Amplitude of the magnetic 
field was varied in interval Bg G [—4,0.95] T with step 
ABq — 0.05 T. The current was increased until the sat- 
uration was achieved. As a criterion of saturation we 
used the relation MjMs > 0.999999, where is the 
total magnetization along z-axis. The resulting depen- 
dence Js{B) in dimensionless form is shown in the Fig. [3] 
by markers. Note a good agreement between theoreti- 
cal prediction and numerical experiment. The reason for 
slight discrepancy in region 6 ^ 1 is that the saturation 
field for the finite-size nanodisk B'^ is slightly smaller than 
Bs. 

It is knowrPSl that the VAL usually appear in pre- 
saturated regime of the ferromagnetic film, see insets of 



FIG. 4: (Color online) Dependence of period of VAL, which appears in the pre-saturated regime, on the applied 
field. The data obtained using micromagnetic modelling (see test) are shown by markers, dark disks corresponds to 
the VAL, whereas the dashed circles represents the fluid-like state of the vortex-antivortex system (see Ref.'^,'^!). The 
solid line corresponds to dependence 2-k/Kq obtained from (18b I. The insets A and B show the magnetization 
distribution in VAL for cases of positive and negative fields respectively. The upper parts of the insets show the 
TO^-component distribution: with decreasing the value of the contour line becomes thinner and when is close 
to the minimum value the contours become dashed. Magnetization distribution within the film plane is shown by 
arrows in the bottom parts of the insets, disks and diamonds show centers of vortices and antivortices respectively. 

The size of the inset area is 50 x 50 nm. 



the Fig. |4] Here we study numerically how the perpen- 
dicular magnetic field changes properties of the VAL. We 
obtained that the positive field (the direction of the field 
coincides with the current direction) increases the con- 
stant of VAL as while the negative field (opposite to the 
current) decreases ag. The resulting dependence ag{b) is 
presented in the Fig. |4] As one can see, the lattice con- 
stant as{b) is very close to the value as = 2i: /Kq, where 
Ko is the wave-vector of unstable magnons for the case 



K'^Xs, see (18b I. Assuming that mismatch between as 



and as remains small for all values of parameters, one can 
use ( 18 1 to obtain the following asymptotical behavior: 
as ~ 1/(1 - 6) for 6 ^ 1 and as ~ l/v1M for h -oo. 



III. CONCLUSIONS 

The perpendicular magnetic field drastically changes 
the process of saturation of magnetic films with spin- 



polarized current. It is shown that the saturation current 
Jg is decreased (increased) in case of codirected (oppo- 
sitely directed) magnetic field and current. There exists 
a critical current Jc > Js which provides "rigid" sat- 
uration - the saturated state which is stable with re- 
spect to the transverse magnetic field of any amplitude 
and direction. The critical current Jc is determined only 
by material parameters (saturation magnetization) and 
thickness o the film. The actions of the perpendicular 
magnetic field and uniaxial anisotropy on the stability 
of saturated state are equivalent. The magnetic field 
changes the constant of the vortex-antivortex lattice as, 
which appears in the pre-saturated regime: as infinitely 
increases if the field approaches the saturation value and 
as decreases if the field is increased in the opposite di- 
rection. For large opposite fields the fiuid-like dynamics 
of the vortex-antivortex system is observed instead of the 
static vortex-antivortex lattice. 



* |alexey@volkov.ca ] 

^ J. Lindner, Superlattices and Microstructures 47, 497 
(2010). 

^ S. Bohlens, B. Kriiger, A. Drews, M. Bolte, G. Meier, and 
D. Pfannkuche, Appl. Phys. Lett. 93, 142508 (pages 3) 
(2008). 

^ A. Drews, B. Kruger, G. Meier, S. Bohlens, L. Bocklage, 
T. Matsuyama, and M. Bolte, Applied Physics Letters 94, 



062504 (pages 3) (2009). 

* A. D. Kent, B. Ozyilmaz, and E. del Barco, Appl. Phys. 
Lett. 84, 3897 (2004). 

5 J. C. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996). 
^ L. Berger, Phys. Rev. B 54, 9353 (1996). 
J. C. Slonczewski, J. Magn. Magn. Mater. 247, 324 (2002). 

* O. M. Volkov, V. P. Kravchuk, D. D. Sheka, and Y. Gai- 
didei, Phys. Rev. B 84, 052404 (2011). 



6 



® Y. Gaididei, O. M. Volkov, V. P. Kravchuk, and D. D. 

Sheka, Phys. Rev. B 86, 144401 (2012). 
" A. Dussaux, B. Georges, J. GroUier, V. Cros, 

A. Khvalkovskiy, A. Fukushima, M. Konoto, H. Kubota, 

K. Yakushiji, S. Yuasa, et al., Nat Commun 1, 1 (2010). 
^ A. Dussaux, A. V. Khvalkovskiy, P. Bortolotti, J. GroUier, 

V. Cros, and A. Pert, Phys. Rev. B 86, 014402 (2012). 



V. Sluka, A. Kakay, A. M. Deac, D. E. Biirgler, R. Her- 
tel, and C. M. Schneider, Journal of Physics D: Applied 
Physics 44, 384002 (2011). 

A. I. Akhiezer, V. G. Bar'yakhtar, and S. V. Peletminskii, 
Spin waves (North-Holland, Amsterdam, 1968). 

^'^ The Object Oriented MicroMagnetic Framework, devel- 
oped by M. J. Donahue and D. Porter mainly, from NIST. 
We used the 3D version of the 1.2Qf4 release, URL http : | 

I //math.nist .gov/oommf /, 



